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We argue that the Standard Model contains stable bound states with a sufficiently large number 
N of heavy quanta - top quarks t,t and gauge bosons W,Z — of the form of collective "bags", with 
a strongly depleted value of the Higgs VEV inside. More specifically, we study one-loop quantum 
corrections to a generic model of them, assuming "quanta" are described by a complex scalar field. 
We follow the practical formalism developed by Farhi et al. for the N = 1 case, i.e. one particle 
in a bag, who found that for a very large Yukawa coupling the classical bags are destabilized by 
quantum effects. We instead study the problem with a coupling constant in the range of the 
Standard Model for a large number of quanta N = 50. .5000. We calculated both classical and 
one-loop effects and found that for such bags quantum corrections are small, of the order of a few 
percent s or less. 



I. INTRODUCTION 

Coherent bound states of a large number of bosonic 
quanta are known as "semiclassical solitons" . An im- 
portant subclass of them - topological solitons - play a 
significant role in physics in general and in the dynamics 
of gauge theories in particular. Electroweak sphalerons, 
for instance, were extensively studied in the context of 
baryon number violation in the Standard Model (SM). 
Instantons drive the breaking of chiral U(l) and SU (Nf) 
symmetries in QCD, while dyons/monopoles are believed 
to be related to the phenomenon of color confinement. 

The objects we will study in this paper belong to the 
less famous subclass, namely "non-topological solitons" , 
which have remained somehow in the shadow of their 
topological relatives, although also studied - as pure the- 
oretical constructions - for a long time. They are lo- 
calized field configurations with finite energy which owe 
their existence not to topological properties, but rather 
to the existence of a conserved charge. The first exam- 
ples were studied by R. Friedberg, T. D. Lee and A. Sirlin 
in [l|, Q and some time later Coleman studied what he 
termed Q-balls 0]. These solitons, as we shall discuss 
below, are built out of scalar fields with an unbroken 
U(l) symmetry. As is well known, Derrick's theorem 
precludes the existence of time-independent soliton con- 
figurations in d > 1 + 1 if only scalar fields are involved. 
However, non-topological solitons, both the Friedberg et 
al. and Q-balls, circumvent Derrick's theorem by an ex- 
plicit time dependence and owe their stability to the con- 
served charge associated to the unbroken U(l). Despite 
of this similarity, they differ in some other respects. For 
instance, in the case of the Friedberg et al. type, there's 
a second symmetry which is spontaneously broken due 
to the presence of a Higgs-like scalar field. They repre- 
sent, as we shall see, regions of space in which the Higgs 
vacuum expectation value (VEV) is suppressed due to 
the presence of heavy particles. In this work we concen- 
trate solely in this kind of Higgs-based solitons and the 



role of quantum corrections as we believe they might be 
relevant in applications to cosmology, although we will 
often compare our results to those previously obtained 
for Q-balls. 

Regarding Q-balls, it was first shown in [§| that they 
exist in the limit of large charge Q, and it was later shown 
in Q that these classical Q-balls actually exist all the way 
down to Q = 1. These small Q-balls may play an impor- 
tant role in cosmological applications [4| and therefore it 
is important to determine the validity of the classical ap- 
proximation. Therefore, the role of quantum corrections 
to these small Q-balls was explored in Q by using the 
methods developed in Q and it was found that although 
they are small in comparison to the classical energy, they 
are comparable to the binding energy for small Q. It 
was finally concluded that quantum corrections (for typ- 
ical values of coupling constants) render small Q-balls 
(Q < 7) unstable. 

In contrast to all these works, we do not discuss hypo- 
thetical objects which can only appear in various exten- 
sions beyond the Standard Model, but address the range 
of stability and existence of the Multi-Quanta Higgs Bags 
(MQHB) in the SM itself. We have in mind the heaviest 
objects of the SM, top quarks t, t and W, Z bosons, to 
which we refer generically as "quanta" . Of course, bags 
with many top quarks, or top -bags as we will call them, 
are quite different from scalar-quanta bags due both to 
their fermionic nature and their coupling to the Higgs 
field. However, these issues, and the potential role of 
these bags in the electroweak baryosynthesis problem, 
will be addressed elsewhere [3, [8( . In this work, which is 
mostly methodical in nature, we first address the issue of 
whether quantum corrections can or cannot destabilize 
these objects. 

The interest in the issue originated from the question 
whether a sufficiently heavy SM-type fermion should ac- 
tually exist as a different state, in which it depletes the 
Higgs VEV around itself and thus is accompanied by its 
own "bag" 0. Although classically this seemed to be 
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possible, it was shown by Farhi et al. [10( that in fact 
quantum (one loop) effects destabilize such bags, except 
at such a large coupling that the theory itself is appar- 
ently sick, with an instability of its ground state. The 
issue rest dormant for some time till Nielsen and Frog- 
gatt [ll| argued that although the Yukawa coupling may 
still be small, there should be binding when the number 
of particles N is such that 
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Indeed, a similar generic argument explains for example 
why gravity, although very weak, can create bound states 
- planets and stars - but only for an "astronomically 
large" N. 

These authors also suggested that 12 top+antitop 
quarks (the "magic number" corresponding to the maxi- 
mal occupancy of the lowest Si/ 2 orbital, 6 quarks and 6 
anti-quarks with all spin and color values) already form 
deeply bound states with near-zero total mass. Here is 
when one of us (ES) became involved: in a paper with 
Kuchiev and Flambaum [12] we checked this suggestion 
and showed that, unfortunately, it is not the case. While 
for a massless Higgs we found a weakly bound state of 12 
tops, for a realistic Higgs mass there is no bound state at 
all. A possible direction to consider is that of introduc- 
ing new hypothetical fermions with a stronger coupling 
(larger masses) than that of the top quark and the possi- 
bility of these forming baryons, well-bound by the Higgs. 
This issue has been discussed further in the second paper 
of the same authors [13j and we will discuss the implica- 
tions of our study for this later. 

Although we disagreed with the calculation by Nielsen 
and Froggatt [ll[ on 12 tops, we think that further stud- 
ies of multi-quanta Higgs bags are well justified. In this 
paper, we will be especially interested in the case in which 
the number of "quanta" N is large enough to consid- 
erably modify the Higgs VEV compared to its vacuum 
value. In other words, we will focus on relativistic bound 
states, in which the "quanta" inside the soliton have 
masses quite different from their vacuum ones, rather 
than non-relativistic atom-like bound states with smaller 
N. This sector of the SM spectroscopy constitutes an 
interesting class of relativistic many-body systems, quite 
different from atoms and nuclei in many respects, espe- 
cially regarding quantum corrections, which we study. 
While the experimental production of a large number 
of heavy quarks/bosons in a sufficiently small volume is 
clearly impossible, at the LHC or any other future pro- 
ton accelerator, such states may play a significant role in 
cosmological baryogenesis at temperatures slightly below 
the electroweak phase transition [8| . 

Turning now to the main subject of this paper; the role 
of quantum corrections to the energy of such configura- 
tions, we express their energy generically as a quantum 
loop diagram expansion 



E(N) = E cla 
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Below we calculate the one-loop coefficient C\ closely fol- 
lowing both the specific model and the practical methods 
developed in Q. Very specifically, the issue we investi- 
gated is the iV-dependence of the 1-loop radiative correc- 
tion. 

The paper is structured as follows. In Section |TT] we 
formulate the Lagrangian and the "thin wall" approxi- 
mation used to prove the existence of solitons at large 
enough N. In Section IIHI we review the methods devel- 
oped in [6t] which provide a practical tool for numerically 
calculating the 1-loop radiative corrections to soliton con- 
figurations. It is important to note that these methods 
were designed for time-independent classical configura- 
tions while the non-topological configurations we just dis- 
cussed have an explicit time dependence by having a non- 
zero N. However, as shown in [5| (in an application to 
Q-balls) the method can easily be adapted to study a 
time-dependent configuration as well. After an estimate 
of the effect of quantum corrections in Section IIVI and 
a discussion of other secondary calculational details in 
Section [V] we proceed to our results in Section fVTl 



II. A CLASSICAL SOLITON 

Consider the action 

(3) 

where <j> is a real scalar field, modeling the SM Higgs, with 
its VEV cf>(r = oo) = v. A complex scalar field ip would 
have "quanta" , our simple version of t, W and Z. This 
system was studied in |l| as an example which exhibits a 
soliton solution at the classical level. In this section we 
review the general arguments which were given in favor 
of the existence of these objects. 

Due to the global symmetry ip -> e la ip in ([3]), there's 
an associated conserved current, d^j^ = 0, with 

j t _ l = i(^*d t ,ip-ipd^ip*) (4) 

and a corresponding conserved particle number 

N = J j°d 3 x. (5) 

Now, let the complex scalar be a time-dependent field 
rotating in internal space with a frequency u>, i.e. 



i/>(x,t) = e- iwt ip(x). 
The classical energy E c i from Q reads, 



(6) 



d 3 3 



(2) 



i(V0) 2 + |v^f+w 2 
A 
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and from ([5]) we have 

N= 2uj J d 3 x\ipf 
The equations of motion read 



(8) 



■ w 2 Vj - v 2 V> + g4> 2 4> = 0, (9) 



■V 2 0-2 5 H 2 0+^(^ 2 - V 2 )</> = 0. (10) 
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By using the equation of motion © for ^ in © and the 
expression (j8]) we find 



E cl = E h + Ncu 



with 



E h 



- Iv/i| 



+ —(h J + 2vhy 

8v z 



(11) 



(12) 



where we have defined as usual h = <p-v. Therefore, the 
classic configuration consists of a static background field 
h(x) and N "quanta" of the -0 field at an energy level 
Eq = w. In the Higgs vacuum, 4>{r) = v, N quanta of 
the ip field would be forced to rest at the bottom of the 
continuum spectrum and the total energy of the system 
would simply be NM. However, in the background of 
a non-trivial Higgs field there are two competing effects. 
On the one hand, the gradient and potential terms in (fT2|) 
increase the energy but, on the other hand, there might 
be some bound states levels with energy < Ej < M 
which can allocate the quanta, lowering the energy of 
the system of particles at the expense of creating such 
a Higgs configuration. In particular, if there's a large 
region of nearly-zero Higgs VEV (where particles would 
be nearly massless) , a soliton will be possible if the energy 
increase by creating such a region is offset by the energy 
decrease of N particles giving away their mass. 

Indeed, the existence of such solitons was shown in [l| 
by considering the following trial configuration: A time- 
independent Higgs field which vanishes inside a sphere 
of radius R while outside it takes its asymptotic value, 
with a transition region of size I. The complex scalar 
ip(t, x) = e~ zut t/j(x) is non-zero inside the sphere and van- 
ishes outside. In a thin-wall approximation, I « R, this 
provides a spherical well in which the complex field forms 
a relativistic bound state, i.e. 



= 0, ip(x,t) = e _iart f simW), r < R, 
(t) = v(l-e-^- R ^ 1 ), r tp = 0, r>R. 



(13) 



with co = ~ . From (fTTjl , (fT2|) the energy of such a config- 
uration, in the limit of large R, reads 



„ nN 7T771 4 o 

E = + R 3 + 0(R 2 ), 

R 2A V h 

which has a minimum at 

3m 4 i? 4 „ 47rm / 3 



(14) 



N: 
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where m 2 = Av 2 /3 and M 2 = gv 2 . Notice the power of 
N with which the energy scales. If the Higgs field took 
its asymptotic value everywhere (i.e. <j){r) = v), the total 
energy of the system would simply be given by E = NM. 
Since E oc TV 3 / 4 « N for large N, there's always a large 
enough N which makes the soliton configuration ener- 
getically favored over simply having N particles in the 
Higgs vacuum. Notice that the critical configuration with 
E c = MN C is given by 



^ IAwm\ 4 
c ~ \ 3M ) 



3_ 
2A' 



(16) 



Since R grows with N, this expression becomes more 
accurate for large N, i.e. when M/m -»■ expression (fl6|) 
gives the limiting value for N c . It was similarly shown in 
[l| that in the limit M/m -> oo, an upper bound on N c 
is given by 



N r < 336 



( m \ 3 1 

\m) A 



(17) 



and therefore when the Higgs becomes massless the soli- 
tons exist not only for large enough N but actually all 
the way down to small N. 

Recall that Derrick's theorem precludes the existence 
of time- independent solitons with an action such as f3j). 
but this no-go theorem is clearly circumvented here be- 
cause of the explicit time dependence in (|13[) . which leads 
to N ~ loA 2 . If we had no time dependence the first term 
in (|14p would vanish, rendering the system's minimum at 
R = 0, in accordance with Derrick's theorem. 

Although we have shown the existence of classical bags 
in certain approximations, in order to explicitly find such 
configurations with minimal energy we adopted a varia- 
tional strategy. We first fix a smooth, 1-parameter, bag- 
shaped, background Higgs field ansatz on which we nu- 
merically solve the equation of motion for ip. We then 
minimize the resulting energy with respect to the varia- 
tional parameter. Once the classical solution is found we 
turn to the study of quantum corrections by implement- 
ing the method developed in Q, which for completeness 
is reviewed in the next section. 



III. QUANTUM CORRECTIONS: 
GENERALITIES 

A. Time-independent configurations 

As we have seen, a non-trivial Higgs background may 
create classical bag-like solitons. However, it also has a 
secondary effect at the quantum level, which is that of 
disrupting the energy levels of the vacuum by the Cas- 
simir effect. To study whether this affects or not the 
existence of solitons we must study the vacuum energy, 
given to leading order in h, by the effective action. The 
calculation of such an object will be plagued by the usual 
divergences of field theory and requires proper renormal- 
ization. In this section we review the method developed 
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by Farhi et al. [6[, which precisely provides us a way 
of calculating the effective action in a manner consistent 
with on-shell mass and coupling renormalization. We 
shall consider the effective action induced only by inte- 
grating out fluctuations in tp. 

The one-loop effective action, induced by integrating 
out the ip field, is given by 



J Dip* D^e lS ^v.a 
and the effective energy is then 



(18) 



E 



eff 



■ 1™ TfSeff, 

1 — >-oo ± 



S ef f = S c i + S c t- logdetff (19) 



where H is the quadratic operator s ^, s ^ , generally of the 

form H = - a +M 2 + V. The matrix identity logdet-ff = 
Tr log H allows us to write the effective action as a trace 
in functional space. We can perform this trace over a 
complete set of functions satisfying Hip a = or 



(-V 2 + M 2 +V)if> a = E a % 



(20) 



which in the free (V = 0) case are simply plane waves 
e lkx . Therefore, in that case, Ef ac ~ / d 4 k\og(k 2 + M 2 ) 
which by standard methods can be brought into the form 

/ d 3 k\/k 2 + M 2 = J d 3 k\E k \, which is the usual vacuum 
energy. In the presence of a background potential there 
might be, in addition to the continuum, some discrete 
levels corresponding to scattering bound states. In such 
case the vacuum energy generally reads 



77 (/ 



f\E a \, 



(21) 



where a = {i, k} labels the discrete (i) eigenvalues, if any, 
and continuous (fc) eigenvalues and E a are the energies of 
the excitations. This is simply the sum of zero-point fluc- 
tuations ^hw where the ^ factor is absent due to ip being 
a complex field. Therefore, there are 3 contributions to 
the total effective energy given by 



E[h] = E cl [h] + E ct [h] + Et ac .[h], 



(22) 



where E ct is the necessary renormalization energy coun- 
terterm. For a fixed h, the spectrum consists of a con- 
tinuous spectrum E 2 = k 2 + M 2 and (possibly) a finite 
number of bound states with energy < E 3 < M , i.e. 

Et ac .[h}= £ E a = Y l E j + T,( 2l + 1 ) [ dkpi{k)E{k) 
a j I J 

(23) 

where pi(k) is the density of states in momentum space 
and the (21 + 1) factor accounts for the degeneracy in the 
angular momentum projection. It can be shown that [ItJ 



it dk 



(24) 



where Si(k) is the scattering phase shift of the Z'th partial 
wave under the potential V(r) and p®(k) is the density of 
states in the absence of potential. Therefore, in order to 
calculate the whole 1-loop effective action we must cal- 
culate the phase shifts Si(k). As is standard in quantum 
mechanics, a Born approximation may be used. This 
represents an expansion in number of insertions of the 
potential V(r) and in this model only diagrams with 1 
and 2 insertions are divergent. Therefore, one defines the 
subtracted phase shift 



^(fc)E5 i (fc)- < 5.( 1 )(fc)-r5, (2) (fc) 



(25) 



where 5 < ^ 1 \k) and (5p^(fc) are the first and second Born 
approximations to Si(k). It is shown in Q that these 
two subtractions can be reintroduced into the action as 
divergent diagrams, which together with E ct [h] combine 
to give a finite contribution denoted by I^f/i]. Therefore, 
we have 

E[h] = E cl [h]+T 2 [h] 

+ i£(2Z + l) / cifc(^(fc))£(fc) + El^l- (26) 

A very useful theorem due to Levinson states that the 
number of bound states ni with angular momentum / in 
a given background potential is given by 



m = -[Si(0)- .Sz(oo)]. 



(27) 



This allows us to partially integrate d/dk in (|26|) . picking 
up a boundary term -M for each bound state, which fi- 
nally gives an effective energy containing 4 terms, namely 

E[h] = E cl [h] + E bound [h] + E cont .[h] + T 2 [h] (28) 

where 



E bomd [h] = -M), 



(29) 



E cont .[h] = --Y,(2l + 1) f^dkStjk)^^ 430) 

7T JO Vfc 2 + M 2 

and E c i is given by (jTTJ) . Notice that after the boundary 
term due to Levinson's theorem is included, the quantum 
correction due to the bound states is negative and there- 
fore it lowers the total energy. This raises the possibility 
that one might find quantum solitons which classically 
are forbidden due to Derrick's theorem. The existence of 
such solitons was explored in Q with no success. From 
the previous discussions, however, we know that classical, 
time-dependent, soliton configurations exist for N > 
and we are interested in determining the effect of quan- 
tum corrections to these classical configurations. 
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B. Time-dependent configurations 

We saw above that the method developed in @ deals 
specifically with time-independent classical configura- 
tions. Therefore, to find the role of quantum corrections 
to the class of non-topological solitons we are interested 
in, we must adjust the methods to a time-dependent con- 
figuration. This has already been done in the case of 
Q-balls in *Sj by one of the same authors who developed 
the methods for the static case. In our case the same rea- 
soning applies and we closely follow the same derivation. 
Nevertheless, as we shall discuss, in applying this method 
to the case at hand there's an important difference with 
respect to Q-balls. The extension to the time-dependent 
case is easily carried out by studying fluctuations in a 
corotating frame. Let 



ip(x,t) 
(j)(x,t) 



(ip c i(x) + r)(x,t)) , 



(pcl(x), 



(31) 
(32) 



with e~ iwt ip c i(x) the classical soliton discussed above and 
Sip = e~ twt ri(x,t) an arbitrary quantum fluctuation. Then 
the classical action reads 



/ 



d A x 



"4! 



(33) 



and to second oder in quantum fluctuations of ip 

Sip = J d 4 x[|(5 -zw)^,t)| 2 -|v77(x,t)| 2 

-g4> 2 d \v(x,t)\ 2 ]. (34) 

Therefore, 

~(d ~iu) 2 + V 2 -g<j> 2 cl (35) 



is the quadratic operator of quantum fluctuations. For 



convenience, we parametrize rj{x,t) 



i(u-E)t 



T)(x), 



so that fluctuations over the static configuration 
{4>ci{x),ip c i{x)} have an energy E - lo and hence their 
contribution to the effective energy is 



Et 



\f a \E a -u\, 



(36) 



where the EaS are the solutions to the eigenvalue prob- 
lem 



{-El 



V 2 + ; 



ll)Va(x) = 0, 



(37) 



i.e. we must study the zero modes of the operator H = 
-E 2 - V 2 + g<p\i- Notice that we are left with a reduced 
problem which is identical to that we would of have to 
solve for the static case (|2TJ)) . with V = g(h 2 + 2hv) = g\ 
and M 2 = gv 2 , but the contribution of each level is shifted 



by lo. Now, since the spectrum in (|37|) is symmetric in 
E -*■ —E, we can write 

Et ac .[h] = ~ Y, \E a - lo\ + \E a + lo\ = Y, max{i; Q , 

^ E a >0 E a >0 

(38) 

As we shall see below, if we only integrate out the fluc- 
tuations Sip from (|3|) we have E a > u> and therefore the 
vacuum energy is again simply 



Etac.[h]= f\E a 

JOL 

and the total energy reads 



E = E cl + £ (\Ei\ - M) + E cont . + T 2 



(39) 



(40) 



with E c i given by (fTTjl . To see that we have E a > lo in this 
case, we note that from the equations of motion we find 
( —lo 2 - V 2 + g<Pli) ipd(x) = in the I = channel, which 
is the most tightly bound state and therefore the lowest 
positive energy is E = \ui\, proving our assertion. It is 
important to note however that had we also considered 
fluctuations of the Higgs field as well, there would be a 
zero mode of the quadratic operator H (which would be 
2x2) now in the I = 1 channel, corresponding to trans- 
lations of the soliton. This would imply that there's an 
energy level Eq below to (in the I = channel), giving an 
additional contribution to the effective energy, namely 
lo - Eq, with respect to (|3"9")l . This is important in the 
case of the Q-ball where one has no other choice than 
considering such a zero mode as there's only one field in- 
volved and actually this term dominates the whole quan- 
tum correction [5J] . This difference is important since this 
additional term, being positive, works against the stabil- 
ity of these solitons and is ultimately the responsible for 
rendering Q-balls with Q < 7, unstable. However, as just 
discussed, that's not the case here due to the absence of 
the translational zero mode in the path integration. 



IV. ESTIMATE OF QUANTUM CORRECTIONS 

Before going to the numerical calculations, we can il- 
lustrate the effects of quantum corrections on the soli- 
tons described in Section [ill by using the thin-wall con- 
struction. This configuration corresponds to a spheri- 
cal potential well of depth -gv 2 = -M 2 and size R, i.e. 
V = -M 2 Q(R - r), and the evaluation of the quantum 
corrections requires solving the energy levels of such a 
system. We wish to find the dependence of these energy 
levels on the size of the soliton. We can easily do this for 
the continuum energy relative to the Higgs vacuum. It is 
given by 

Econt. = J d 3 x J d 3 k [Vfc 2 + M 2 + V-Vk 2 + A/ 2 ] 



(41) 
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and therefore 



Econt. = (4tt) 2 J r 2 dr J dkk 2 [fc - Vfc 2 + M 2 ] . 

The integral over k is obviously divergent and requires 
proper regularization. However, the ^-dependence is 
given by 



(42) 



and therefore it grows as the volume of the soliton and 
its effect is that of modifying the coefficient of TV 3 / 4 in 
the classical energy (|T5|) . The precise dependence of the 
bound states' contribution is harder to find analytically 
and will depend strongly on the coupling g, as this deter- 
mines the depth of the potential -gv 2 . Nevertheless, we 
can understand its main effect. We know that a spher- 
ical well of depth -gv 2 has a minimum size Ro = iv^jg 
at which the first bound state appears. For large cou- 
pling our potential becomes deeper and a large number 
of states could is bound. As is well known, in the case 
of an infinite spherical well there are an infinite number 
of bound states given by the spherical Bessel functions 
ji(kr). The boundary condition ji(kR) = determines 
therefore the energy levels of a relativistic particle to be 



E 



ill 



1 

R 



(43) 



where (3 n i represents the n'th zero of the Vth order Bessel 
function. The number n(R) of bound states in a finite 
potential well is finite and grows with R and therefore 
the contribution ([29]) will be given by n(R)(l/R - M), 
which may be significant at large R. Indeed, as we shall 
comment more later, for large enough coupling, it may 
become so negative that it wins over the classical energy 
of the soliton, making the vacuum unstable under the 
creation of empty bags of zero Higgs VEV. For relatively 
small values of g (~ 1 - 3) this term has the finite effect 
of lowering the number N c of particles required to have 
a stable configuration. 



V. CALCULATIONAL DETAILS 
A. Variational method for classical solution: 

Quantum corrections require some operations with the 
background Higgs field (such as finding its Fourier trans- 
form) which is difficult in practice to do for numerically 
generated shapes. Therefore we adopted a variational 
approach, in which the field is given in a simple analytic 
form, with certain parameters to be determined from en- 
ergy minimization. 

Our purpose is to find a variational minimum of (|40|) . 
For this we used a family of 1-parameter Gaussian back- 
grounds for the Higgs 



ft(r) 



-ve 



-r 2 v 2 /2i 



(44) 



where w > is the variational parameter representing the 
size of the configuration. Notice that ((f)) ~ for r < w 
while (4>) ~ v for r > w. As is customary, we refer to this 
configuration as a bag. It represents a region of size ~ w 
of symmetric phase embedded in a space of broken phase. 
Inserting this ansatz into (|12|) gives us an analytical ex- 
pression for the Higgs energy, which we denote Eb ag (w). 
Notice also that this profile for the Higgs provides an at- 
tractive potential gx( r ) under which the field if) scatters. 
By solving equation (|37[) with I = we find ui, which is 
the lowest bound state level, henceforth denoted by E . 
Recall that the energy of the soliton is given by 



E(w) = E bag (w) + NE (w). 



(45) 



For the energy to have a minimum we require that 
(SE)n = and therefore 



E>(w)+NE' (w) = 



(46) 



with ' denoting derivative w.r.t. w. This defines A'* as a 
function of the bag's size w by 



E' a (w) ' 



(47) 



Furthermore, we must also have E bag (w) + NE (w) > 0. 
By taking the derivative of the first equation and noting 



that E n < 0, we find 



N'(w) > 



(48) 



to have a minimum. We are therefore interested to see 
if there's a range for which N'(w) > with N given by 
63). 



B. Quantum Corrections 



As explained in the Introduction we will follow here 
[(| , who developed a rather practical way to calculate the 
quantum correction, due to one-loop diagrams involving 
the "matter field" ip, to the total energy of the system. 
It is given by several terms; the binding energy of the 
bound states, -Efeound, the continuum term, E cont ., repre- 
senting regularized scattering phases, and the remaining 
regularized mass and charge renormalization, T2- 

To evaluate Etound we proceed as follows: Recall that 
Si(oo) = and therefore, according to Levinson's theorem 
(|27p , the number of bound states with angular momen- 
tum I is given by 



1 



nj = -*,(* = 0). 



(49) 



According to @, from the radial equation of motion, 
one finds that the phase shifts are given by S[(k) = 
-2mPi(k,r = 0) with the #'s satisfying 



- 2ikpi(kr) + 2(A) 2 + -g X (r) = (50) 
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where pi(x) = d/dx (in [x/i^x)]) and hj(x) is a spherical 
Hankel function and ' denotes derivative with respect to 
r. By solving this differential equation and taking the 
limit k -*■ we found the number of bound states to be 
expected in every angular momentum channel at a given 
w. With this knowledge, we used a shooting method to 
find the energies of such states by solving the differential 
equation 

(-V 2 +g^ cl )^ = E 2 ^. (51) 

That is, we repeatedly solved this differential equation 
for ip = — in a box of size L = 30 for increasing values 
of E from to M, with boundary conditions u(0) = 0, 
it'(0) = 1. We then selected those solutions which vanish 
at infinity (i.e. u(L) = 0), verifying that the number 
of such solutions matched those predicted by Levinson's 
theorem. 

To evaluate the continuum contribution (|30[) , one must 
find the Born approximation to the phase shifts. These 
are similarly found by solving a set of coupled differential 
equations, which are detailed in the Appendix. 

The calculation of the renormalization term T2 is de- 
tailed in the Appendix and involves evaluating a dou- 
ble integral which is readily done numerically. It will be 
shown however that it is numerically suppressed to be 
several orders of magnitude smaller than the other con- 
tributions and may be neglected. 



VI. RESULTS 

We discuss the results of our numerical calculations 
starting with the lowest s-wave bound state. Its depen- 
dence on the bag size w is shown in Fig. [T] The points 
represent the numerical results and the continuous line 
is an interpolation which, from (|47p , allowed us to find a 
smooth expression to be used below. Note that the curve 
does not go through the w = 0.9 point, in which case there 
is no bound state: the point simply shows the mass of the 
quantum outside the bag, E = M = 1 for the coupling 
we use. Obviously, there's a minimum bag size u> m in ~ 1 
at which the first state gets detached from the continuum 
to become a bound state with an energy just slightly be- 
low M. It is also clear that such shallow levels should 
be described by an atomic-like approach rather than by 
collective bags: we will not investigate their properties 
in this work, focusing only on deeper bound ones. 

Table Q] shows some results for the energy of the bound 
states in different /-channels for different sizes of the bag. 
As we increase w, new levels appear and the energy of the 
lowest bound state approaches 0. We also observe that 
the levels appear in the same order as in the spherical 
well potential (see, e.g. [111). For instance, for w = 7.9 
the order of the energy levels is Is, lp, Id, 2s, 2p, 2d, 3s. 
The number of quanta to bind in order to ensure that a 
bag of width w has a saddle point is shown in Fig. [5] It 
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TABLE I: This table shows a typical set of bound state levels 
in different angular momentum channels for different sizes of 
the bag with g = A = 1. 
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w 

FIG. 1: Energy of the lowest bound state as a function of 
w for g = X = 1. The points show the results of our numeri- 
cal calculations while the continuous line is an interpolation 
which allows us to find a smooth expression for TV. We can 
see that there's a minimal value w m i n ~ 1 at which a bound 
state appears and its energy decreases smoothly to zero as 
w increases. The number of levels in each channel coincides 
with those predicted by Levinson's theorem. 



continuously increases with w: as we showed in the previ- 
ous section the condition N'(w) > needs to be satisfied 
in order for the saddle point to be a local minimum and 
as we can see from the figure, this is always the case for 
w > w m i n . By taking g = A = 1 in equation (|16p . we ex- 
pect a stable soliton to be found for N ~ 51. Indeed, our 
results in Fig. [3] show that without taking into account 
quantum corrections N c ~ 52 and the effect of quantum 
corrections is to lower N c to ~ 50. Recall that we are 
ignoring the effect of the Higgs self-coupling at the quan- 
tum level and therefore A enters only in the expression 
for the classical energy. By taking m/M -»■ (i.e. a BPS- 
like limit), the Higgs becomes massless and therefore the 
exchange forces between quanta become Coulombic. In 
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FIG. 2: Number of particles necessary to extremize the energy 
of a configuration of bag's size w. The number of particles 
to create a minimal bag is ~ 20, with a size w m in ~ 1. The 
dashed line represents the asymptotic behavior ~ w 2 ' 86 for 
large w. 



FIG. 3: This figure shows the binding energy of the soliton 
with g = 1 as a function of the number of particles. The solid 
line represents the binding energy including the quantum cor- 
rection, while the dashed line contains only the classical en- 
ergy. We can see that after a critical point N c - 52, it is 
energetically favorable to bind particles in a finite size bag 
rather than contracting the bag to zero size forcing the parti- 
cles to migrate to the lowest state of the continuum spectrum 
with energy JVM. 



this limit the minimum number of quanta N c to form a 
bound state is expected to be significantly decreased, as 
confirmed by equation (fT?)) . Notice also that this limit 
can be taken while keeping g, and therefore quantum cor- 
rections, small. Our numerical results for A = 10 and 
g = 3 are shown in Fig. 3] By further increasing the cou- 
pling g while keeping A ^ the size of the critical bag can 
be significantly decreased. For what we have seen, quan- 
tum corrections would not seem to spoil the existence 
of such small solitons but rather promote it, contrary to 
the case of Q-balls, due to the absence of the translational 
zero mode as discussed in §111 Bl 

The vacuum energy, relative to the classical energy is 
embodied in the coefficient C±(N), which is shown in 
Fig. [5] Based on our discussion in Section [TV] we expect 
Ci(N) = |f ^° c - to behave as a - for large N, with 
a determined by the (renormalized) contribution from 
the continuum, while b and f(N) are determined by the 
bound-states contribution. Indeed, we found that our 
results up to N ~ 6000 are nicely fit by a = 1.64, b = 2.04 
and a very small power f(N) = TV 03 : A logarithmic fit 
is also possible. Since this plot is for g 2 /47r « .08, we 
conclude that the one-loop coefficient G\ in the generic 
expansion ([2]) grows slowly with TV and is asymptotically 
of order 1. Therefore, even for the rather large coupling 
of top quarks (g t ~ 1), quantum corrections seem to be 
under good theoretical control. 



VII. SUMMARY AND DISCUSSION 

The main content of the present paper is the explicit 
calculation of the regularized one-loop quantum correc- 
tion to the energy of Multi-Quanta Higgs Bags, in a 
model in which "quanta" are represented by a charged 
scalar field. The main results are: 

(i) For the range of couplings corresponding to the heav- 
iest fermion (the top quark) and gauge bosons of the SM 
such bound states exist for N larger than about 50; 

(ii) quantum corrections to their energy is reasonably 
small, at the level of few percents at large N, although 
they play a larger role at smaller N, where they tend to 
lower the binding boundary N c . We do not think that 
the region close to it is under firm theoretical control, 
unlike the domain of large N > 100. 

(iii) For large solitons, the main quantum correction is 
given by the scattering states continuum, which grows 
with N as the classical energy itself. 

(iv) Even extrapolated to very large N we find that the 
limits of validity of the perturbative expansion do include 
g « 1, i.e. the top quark of the SM. 

As mentioned in the Introduction, the modification of 
the (quantum bound states and classical bag) calculation 
for fermions (real top quarks) and gauge bosons W,Z of 
the SM will be carried out elsewhere together with 
cosmological applications [1] . 

As a remaining discussion issue, we may end up with 
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FIG. 4: This figure shows the binding energy of the soliton 
with A = 1CT 4 and g = 3 as a function of the number of par- 
ticles. Again, the dashed line is the classical binding energy 
while the solid line includes the quantum corrections. The 
number of particles to stabilize the bag is decreased by a 
smaller Higgs self-coupling and quantum corrections still fa- 
vor the existence of such a soliton. 
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FIG. 6: All 4 contributions to the total energy in a log scale. 
Ebound and F2 are negative so their absolute values are shown. 
We can see that the main quantum correction for small soli- 
tons is given by Ebound while for large solitons E cont . becomes 
the most relevant. The renormalization term T2 is numeri- 
cally suppressed by several orders of magnitude. 




2000 4000 6000 8000 10000 12000 14000 

N 



FIG. 5: Coefficient Ci(N) for g = A = 1. The dashed line 
represents an extrapolation of our results (solid line). 



the question when and what happens if the coupling is 
larger than that of the top quark, namely g > 1. It has 
been argued by Farhi et al. [6| that for sufficiently large 
g the model becomes sick, because the vacuum itself (i.e. 
the N = state) becomes unstable against the sponta- 
neous creation of large empty bags. Our experience with 
the large-coupling domain confirms it. A closer look at 



each contribution to the vacuum energy for larger g un- 
covers the reason for this instability. From Figs. [5] and [7] 
we can discern the contribution to the total energy from 
each separate term in the effective energy for g = 1 and 
<7 = 3, respectively. For g = 1 we observe that T2 is numer- 
ically suppressed for all N and the main quantum cor- 
rection is due to an interplay between E uoun d and E cont ., 
the former being dominant (and negative) for small N 
(< 100), while the latter becomes dominant (and positive) 
for large N and behaves similarly to the classical energy, 
as expected from our discussion in Section IIVI Already 
for g = 3 we found a significant increase in the contri- 
bution E uoun d, which becomes comparable in magnitude 
(but opposite in sign) to the classical energy. Note that 
such coupling would correspond in the SM to fermions 
with a mass ~ ,5TeV , and e.g. baryons made of such hy- 
pothetical fermions were suggested in [l3j |. We conclude 
that those states, if such fermions do exist, would have 
noticeable quantum corrections. Thus their theoretical 
status is much less solid than that of top-balls, perhaps 
asking for dedicated further study. 

For g > 5 the bound energy term under discussion be- 
comes so negative that it completely compensates the 
classical energy of even an empty bag, therefore render- 
ing the vacuum unstable. However, at such strong cou- 
pling the theory is in serious trouble and clearly beyond 
the range of validity of any perturbative expansion. 

Acknowledgments. 
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FIG. 7: All 4 contributions to the total energy for g = 3. For 
larger coupling we can see how the contribution from Ebound 
increases roughly an order of magnitude. It becomes com- 
parable in magnitude to E con t. (although with the opposite 
sign) for large N and dominant at small N. 



Appendix 



Here we summarize the equations needed for the cal- 
culation of the Born approximation to the phase shifts. 
This is needed for the explicit calculation of the contin- 
uum contribution. We denote the first and second Born 
terms by Sn and A series expansion in powers of g 
for the phase shifts may be written as 

p l =gf3 ll +g 2 [3 l2 + ... (52) 

with 6u = -2gRefin(k,r = 0) and Si 2 = -2g 2 Re(3i 2 (k,r = 
0) and /3;i,2 satisfy the set of coupled differential equa- 
tions 

-ifi 1 -2ikp l (kr)p' ll + ±x(r) = 0, (53a) 
-2ikpi(kr)P l2 + 2(p' ll ) 2 = 0. (53b) 



>12 



The continuum contribution is given by 

e cont .(k)dk, (54) 
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FIG. 8: Typical fc-dependence of the Born-subtracted contin- 
uum energy for 3 values of increasing w from the bottom up. 
The integral over k grows quickly with increasing w. 



with 



Scontik) = ~- £(2Z + l)fr(fc) (•*) 



Vfc 2 + M 2 



and Si(k) = 5i{k)-5} '{k)-5} (k). Each of these terms is 
given as a solution to the set of differential equations (foTJj) 
and (|53[) . In figure ([5]) we show the continuum energy 
£cont. as a function of k for 3 increasing values of w from 
the bottom up. It can be seen that the Born subtracted 
61 is rapidly convergent in k and that the integral grows 
quickly with w. 

The renormalization effect, encoded in P2 is calculated 
from a set of Feynman diagrams and is given by 



J 2tt 2 H Jo 



log 



(4tt) 2 
M 2 + q 2 x(l- x) 



dx 



2„2 



4x(l - x)v q 
M 2 - m 2 x(l - x) 



h{qf 



m 2 x(l - x) 



M 2 - m 2 x(l - x) M 2 - m 2 x(l - x) 



where x(l) an d h(q) are the spatial Fourier transforms of 
x(x) and h(x). Notice that there are numerical factors 
(47t) 4 in the denominator, which suppresses this contri- 
bution with respect to the other ones. See [6j for more 
details. 
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